Re normalisation and fixed points in Hilbert Space 
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The energies of low-lying bound states of a microscopic quantum many-body system of particles 
can be worked out in a reduced Hilbert space. We present here and test a specific non-perturbative 
truncation procedure. We also show that real exceptional points which may be present in the 
spectrum can be identified as fixed points of coupling constants in the truncation procedure. 
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The ab initio construction of a rigorous quantum 
many-body theory able to describe bound particle sys- 
tems like molecules, atoms, aggregates, atomic nuclei and 
condensed systems has developed over a long period of 
time starting in the sixties Q. 

In practice the explicit resolution of the problem ne- 
cessitates the diagonalisation of the many-body Hamilto- 
nian in Hilbert space which is spanned by a complete set 
of basis states, in principle of infinite dimension, at least 
generally very large. In many cases the information of 
interest is however restricted to the knowledge of a few 
energetically low-lying states which often possess collec- 
tive properties. Hence it would be convenient to work in 
a finite truncated subspace of the original Hilbert space. 
Many different procedures have been proposed and suc- 
cessfully applied in the framework of the nuclear shell 
model and related microscopic descriptions. 
Rigorous projection methods lead to effective Hamilto- 
nians which can in principle be explicitly generated by 
means of perturbation techniques 0, 0] ■ Unfortunately 
there exists no straightforward control on the conver- 
gence properties of the perturbation expansions which 
are involved, especially not when the interaction between 
the particles is strong as it is the case in atomic nu- 
clei or in quantum spin systems for instance. Many at- 
tempts have been made in order to overcome this prob- 
lem |3j, |4j, |5|, |6( . Pragmatic phenomenological procedures 
have also been introduced 0, 0- More recently effec- 
tive 2-body interactions have been constructed by means 
of a non-perturbative renormalisation technique which 
cuts-off the large momentum components of the interac- 
tion 0. 

The investigations which follow are developed in the 
spirit of former work based on renormalization con- 
cepts ^1 Ell 10 ■ We introduce here a possible 
non-perturbative scheme for the bound state many-body 
problem which relies on a systematic reduction of Hilbert 
space. 

Formal framework. Consider a system with a 
fixed but arbitrary number of bound quantum ob- 
jects in a Hilbert space TiS^ of dimension N gov- 
erned by a Hamiltonian i?W \ dy*\ 9%^^ 5 ■■■9p N j where 



{9%^ 1 92*\ ' " > 9p^ *—> 9^} are a set of parameters 
(coupling constants) which characterise ff'*' . The eigen- 
vectors \M N) (g m )) {i=l,---,N} of tfW span the 
Hilbert space and are the solutions of the Schrodinger 
equation 

ffW (5 W)l*w (5 W)) = . (i) 

The diagonalisation of IfW delivers both the eigen- 
values {Aj^W), i = 1, • • • ,N} and eigenvectors 

i = 1, • • ■ ,n\ in terms of a linear com- 
bination of orthogonal basis states i = 1, • • • , N}. 
Since dim = N is generally very large if not infinite 
and the information needed reduces to a finite part of 
the spectrum it makes sense to try to restrict the space 
dimensions. If the relevant quantities of interest are for 
instance M eigenvalues out of the set {^i(g^ N ')\ then 
one looks for a new Hamiltonian H( M \g( M }) such that 

= W A '))|*fV A '))> (2) 

with the constraints 

U9 iM) )=U9 {N) ) (3) 

for i = 1, ,..,M, Eq. (3) implies relations between the 
sets of coupling constants g( M ^ and g^ N ^ 

gk - Jk[9i i92 > -9p ') ( 4 ) 

with k = I,..., p. The solution of this set of equations 
generates the effective Hamiltonian H^ M ' whose spec- 
trum in the restricted space H (M ^ is the same as the 
corresponding one in the original space. 

Renormalisation algorithm for systems at temperature 
T = 0. We develop an explicit but general approach 
which allows to implement the former procedure by fol- 
lowing the evolution of the effective Hamiltonian of the 
system when the dimensions of the Hilbert space are 
systematically reduced. Using the Feshbach formalism 
[15J, |l3 we divide the Hilbert space 7v- N > into two sub- 
spaces, PH {N) and QH {N) with 

dim PH {N) =N-1, dim QH {N) = 1 . (5) 
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In the projected subspace PTi^ the system with energy 
E is described by the effective Hamiltonian 

H eff (E) = PHP + PHQ(E - QHQy 1 QHP . (6) 

H is written in the general form 

H = Hq + gH x (7) 

where Hq and Hi are Hamiltonian operators and g a 
strength parameter (coupling constant) which takes the 
value gW in «W. 

The complete set of basis states i = 1, ■ • • , iV} 

may f.i. be chosen as the eigenvectors of Hq with the cor- 
responding eigenvalues {e^, i — 1, • • • , TV}. The expres- 
sion iJ e / / (E) is generally the starting point of theories 
which use perturbation expansions |l| . Here we proceed 
differently. We consider 

N-l 

i=l 

which is the projection on PH^' of an eigenvector 

N 

iM JV) )=E a S°(ff w )i*«> (9) 

i=l 

of HS N \ If is the eigenvalue corresponding to |\E r ] iV ^) 
we look for the solution of 

F e// (Af >)P|*J">) = AWP|*W> . (10) 

We consider P|<I'^ V - ) ) to be the lowest energy eigenstate 
and QW N ' to contain f. i. the highest one in energy. 
Any other state may in principle be chosen. We impose 
the lowest eigenvalue in the PHS N } subspace to be the 
same as the one in the complete space 

Projecting the expression of Eq. (|10|) on ($! | which is the 
eigenvector of H with lowest energy 

($^ e //(Ar ) )ip*w ) =Ar ) ( ff w )4f ) ( 5 (Ar) ) . (12) 

and defining the Hamiltonian which acts in the projected 
space PH {N) as H^ N ~^ = H +g^ N -^ H u Eq. JTl} leads 
to 

(^\H eff (x{^)\P^)=^g(^). (13) 

where 

S[g K ')-H 1N +H 1N (A 1 -H NN ) H N1 

(14) 

with 

= (15) 



and 

H^={<P 1 \H^\P^[ N) ). (16) 

H^[ is the same expression as H^ ^ with ($i| re- 
placed by ($at|. 

The r.h.s. of Eq. (|13|l can be worked out explic- 
itly. The denominator in the second term of H e f f (A^ ) 
of Eq. (|14|) is a scalar quantity since dimQWW = 1 
where N designates the highest-lying energy eigenstate 
of if , |$at). Imposing the constraint introduced through 
Eq. Hll|l leads to a relation which fixes g^ N ^', One gets 
explicitly 

a tN-l) g {N-X)* + b (N-l) g (N-l) + C (N-1) = 

where 

a (jv-i) = Gin _ HnnFin (i 8 ) 

with 

^■ = (^1^1$,) (19) 

= a{VH NN (\{ N) ~ ex) + F 1N (X[ N ^ e N ) (20) 

c (iV-i) = _ a W (A W_ ei)(A W_ ejv)) (21) 

with 

N-l 

F 1N =J2 a H ) (*iWi\$i) ■ (22) 

i=l 

and 

N-l 

Gin = H 1N ]T «if } (*jv|^i |^) . (23) 

i=l 

Since Eq. I|17fl is non-linear in g^^ 1 ^ and has two so- 
lutions, g( N ~^ is chosen as the one closest to g( N ' by 
continuity. The process can be iterated step by step by 
projection from the space of dimension N — 1 to N — 2 
and further, generating subsequently a succession of val- 
ues of the strength parameter (coupling constant) g( k > at 
each iteration. At each step the projected wavefunction 
\P^^) is obtained from ) by elimination of a state 
|$fe). Taking the continuum limit for large N leads to a 
differential flow equation for g. 
Example. 

We consider a real symmetric tight-binding Hamilto- 
nian with Ho = 0, diagonal elements = fj and 
non-diagonal ones ($i\Hi\$ i+ i) = = 7 
which is typical of strongly coupled systems. By essence, 
a reduction of Hilbert space by means of perturbation 
expansions is meaningless in this context. Starting 
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with an initial Hilbert space dimension N we apply 
the renormalisation procedure described above to g 
starting from an initial value jW, The evolution of the 
lowest eigenvalues and the flow of g are shown in Table 
1 for different values of g and N. Several conclusions 
can be drawn. First, as expected, the stability of the 
spectrum of low-lying states is the better the smaller 
the non-diagonal coupling between neighbouring states. 
Second, the stability increases when the initial dimension 
of the space increases. Third, the coupling constant may 
considerably change, decreasing in the present case. As 
a general prescription for the algorithm to work one 
should always keep the states which are strongly coupled 
those basis states which contribute in an essential way 
to the low-lying states of the physical spectrum. 



N 


n 


g 


Ai 


A 2 


A 3 


A 4 


A 5 


10 


10 


20 


0.81 


3.17 


6.90 


11.69 


17.15 




7 


13.4 


1.02 


3.93 


8.29 


13.43 


18.57 




5 


8.18 


1.10 


4.09 


8.18 


12.27 


15.26 


20 


20 


20 


0.22 


0.89 


1.98 


3.47 


5.34 




10 


3.28 


0.26 


1.04 


2.26 


3.83 


5.62 




5 


1.13 


0.30 


1.13 


2.25 


3.38 


4.21 


30 


30 


20 


0.10 


0.41 


0.92 


1.62 


2.51 




15 


6.02 


0.12 


0.46 


1.02 


1.77 


2.68 




5 


1.03 


0.14 


0.52 


1.03 


1.55 


1.93 


50 


50 


20 


0.04 


0.15 


0.34 


0.60 


0.94 




20 


3.73 


0.04 


0.17 


0.37 


0.65 


0.99 




5 


0.38 


0.05 


0.19 


0.38 


0.57 


0.71 


20 


20 


1.0 


0.01 


0.04 


0.10 


0.17 


0.27 




10 


0.33 


0.01 


0.05 


0.11 


0.19 


0.28 




5 


0.11 


0.015 


0.06 


0.11 


0.17 


0.21 



Table 1. 

Evolution of the coupling constant and the 5 lowest 
eigenvalues of the tight-binding matrix described in 
the text. Here = 1, 7 = 0.5. iV is the initial space 
dimension, n the restricted dimension and g the running 
coupling constant. 

Exceptional points and fixed points. 
In the present renormalisation scheme exceptional points 
which generate divergences in perturbation expansions 
are related to the existence of fixed points of the coupling 
constant g. 

It has been rigorously established that the eigenvalues 
Afc(<?) of H{g) = H + gHi are analytic functions of g 
with only algebraic singularities 0, 0, 0] . They get 
singular at so called exceptional points g = g e which are 
first order branch points in the complex g - plane. Branch 
points appear if two (or more) eigenvalues get degenerate. 



This can happen if g takes values such that Hkk = Hu 
where Hkk = (&k\H\&k) which corresponds to a so-called 
level crossing. As a consequence, if a level belonging to 
the PTL subspace defined above crosses a level lying in 
the complementary QH subspace the perturbation devel- 
opment constructed from H e ff(E) diverges 11(1 . Excep- 
tional points are defined as the solutions of |l8| 



f(X(g e )) = det[H(g e ) - X(g e )I] = 



and 



#(%)), 
dX 



!a=a( Sc ) 







(24) 



(25) 



where f(X(g)) is the secular determinant. It is now pos- 
sible to show that exceptional points are connected to 
fixed points corresponding to dg/dx = in specific cases. 
If {Xi(g)} are the set of eigenvalues the secular equation 
can be written as 



N 



n (a - a,-) = 



(26) 



i=i 



Consider X — X p which satisfies Eq. (j21J). Then Eq. f2"5|) 
can only be satisfied if there exists another eigenvalue 
A<? = A p , hence if a degeneracy appears in the spectrum. 
This is the case at an exceptional point. 

Going back to the algorithm described above consider 
the case where the fixed eigenvalue Ai gets degenerate 
with some other eigenvalue x\ k \g = g e ) at some step k 
in the space reduction process. Since Ai is constrained 
to be constant, 



Aj (fe) (5e) = ^'(g'e) 



(I), 



(27) 



which is realised in any projected subspace of size k and 
I containing states |4>i) and |<I>i). Going over to the con- 
tinuum limit for large values of N and considering the 
subspaces of dimension x and x + dx in this limit one can 
write 



dXi 
dx 



= = 



dXj(x) 
dx 



Consequently 



dXj dg e 
dg e dx 







(28) 



(29) 



Due to the Wigner - Neumann avoided crossing rule the 
degeneracy of eigenvalues is generally not fulfilled for real 
values of the coupling constant and the derivative of A^ 
with respect to g vanishes. There exist however specific 
situations, like systems with special symmetry properties 
[ToL l2~f| or infinite systems 01 f° r which degeneracy for 
real g can occur. In these cases Eq. (|29[l is realised if 



dge 
dx 



and 



dX t 
dge 



^0 



(30) 
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The second relation works if crossing takes place and g e 
is a fixed point in the sense of renormalisation theory. 

Eq. H^Ojl shows the connection between exceptional 
and fixed points in the framework of the present ap- 
proach. Ground state degeneracy due to level crossing 
is indeed a signature for the existence of phase transi- 
tions [20|, perturbation expansions break down at these 
points. The ground state wavefunction changes its prop- 
erties when the (real) coupling constant g crosses the 
exceptional point g e . There the eigenstates exchange the 
main components of their projection on the set of ba- 
sis states = 1, ...N. It is worthwhile to notice that 
the same result is not restricted to the ground state, it 
is valid at any physical level crossing in which one of the 
eigenvalues stays constant for any value of the coupling 
constant. We have tested this property on several sys- 
tems. 

Example. 

Consider an N x N matrix which corresponds to a 
highly degenerate and strongly coupled system located 
at a fixed point whose representative matrix possesses 
equal diagonal matrix elements {Hq = 0, Ha — —0.5 in 
the present example) and equal non-diagonal elements 
(Hij — +0.5). Hij is defined in Eq. l(T?ty. The diagonali- 
sation and the application of the Hilbert space reduction 
procedure leads to rigorously stable lowest eigenstates 
equal to —20 and the coupling constant which charac- 
terises the system stays at its initial value, here g = 20, 
for any value of N. The system is located at a fixed 
point corresponding to a state level crossing (exceptional 
point) as discussed above. 

Conclusions. 

In summary, we developed a non-perturbative effective 
theory of the bound state many-body quantum prob- 
lem based on a reduction process of the dimensions of 
the initial Hilbert space. The central point concerns the 
renormalisation of the coupling constant which charac- 
terises the initial Hamiltonian under the constraint that 
the lowest eigenenergy which corresponds to the ground 
state of the system should not change. We showed the 
relationship which exists between exceptional points cor- 
responding to level crossings in the spectrum where per- 
turbation expansions break down and fixed points of the 
coupling constants which characterise phase transitions. 
The approach was tested and discussed on explicit exam- 
ples of strongly coupled and highly degenerate systems 
encountered in condensed matter physics. The formalism 
can be applied to other physical quantum systems like 



atoms, molecules, aggregates as well. It can be extended 
to systems characterised by several coupling constants 
and at finite temperature. Effective operators acting in 
reduced space can be worked out. We shall present these 
developments in forthcoming work. 
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